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ABSTRACT 



The multipole moments of the power spectrum of large scale structure, observed in redshift 
space, are calculated for a finite sample volume including the effects of both the linear velocity 
field and geometry. A variance calculation is also performed including the effects of shot noise. 
The sensitivity with which a survey with the depth and geometry of the Sloan Digital Sky Survey 
(SDSS) can measure cosmological parameters fio and 60 (the bias) or Ao (the cosmological 
constant) and bo is derived through fitting power spectrum moments to the large scale structure 
in the linear regime in a way which is independent of the evolution of the galaxy number density. 
A fiducial model is assumed and the region of parameter space which can then be excluded to a 
given confidence limit is determined. In the absence of geometric and evolutionary effects, the 
ratios of multipole moments (in particular the zeroth and second), are degenerate for models 
of constant f3 w i7°-^/6o. However, this degeneracy is broken by the Hubble expansion, so that 
in principle Qq and bo may be measured separately by a deep enough galaxy redshift survey 



(Nakamura, Matsubara, & Suto (1997)). We find that for surveys of the approximate depth 
of the SDSS no restrictions can be placed on JIq at the 99% confidence limit when a fiducial 
open, flo = 0.3 model is assumed and bias is unconstrained. At the 95% limit, ilo < -85 is 
ruled out. Furthermore, for this fiducial model, both fiat (cosmological constant) and open 
models are expected to reasonably fit the data. For flat, cosmological constant models with a 
fiducial r^o = 0.3, we find that models with ilo > 0.48 are ruled out at the 95% confidence limit 
regardless of the choice of the bias parameter, and open models cannot fit the data even at the 
99% confidence limit. We also find significant deviations in f3 from the naive estimate for both 
fiducial models. Thus, we conclude for the SDSS that linear evolution- free statistics alone can 
strongly distinguish between Clo — 1 and low matter density models only in the case of the 
fiducial cosmological constant model. For the open model, JIq = 1 is only at best only nominally 
excluded unless < 0.3. 



Subject headings: cosmology: Large Scale Structure of the Universe 
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1. Introduction 

With the expectation of new and much larger samples of galaxies with measured redshifts, e.g. the 



Sloan Digital Sky Survey (SDSS, see e.g. Gunn & Weinberg (1995), Strauss (1997)), it is opportune to 



reconsider what information may be extracted from these data. In particular, we would like to determine 
how well fundamental cosmological parameters like the current matter density f^o a-nd cosmological constant 
Aq can be inferred from redshift surveys. Perhaps the simplest approach would be to just measure the 
mean number density of galaxies as a function of redshift. Unfortunately, this number density is evolving 
through processes other than cosmological expansion, and these factors must be accurately deconvolved if 
we are to extract the desired information about cosmology. Alternatively, in the standard cosmological 
model, the growth of perturbations in galaxy number density is driven only by gravity, so we expect that 
measurements derived from these variations will be cleaner signals of cosmology. Locally, this perturbation 
growth induces peculiar-velocities in the galaxies which lead to distortions in redshift maps. At increasing 
redshift, the effects of these distortions evolve due to the changing growth rate of perturbations, and due to 
the geometry and expansion of the universe. This evolution, in turn, depends on fundamental cosmological 
parameters like the curvature, fioj f^nd Aq, so one may hope to infer values of these parameters by observing 



the change in the distortions with redshift. Nakamura, Matsubara, & Suto (1997) have recently calculated 



the effects on the expectation value of the various multipoles of the correlation function, to first order in the 
redshift. However, they made no estimate of the expected errors, necessary for determining the sensitivity 
of any real survey. 

In this paper we shall derive the mean, cosmic variance and shot noise associated with the statistics of 
multipole moments of the linear power spectrum calculating the redshift-dependent effects to all orders. 
Note that moments beyond the zeroth are induced by redshift distortions and would not be present in 
the mean for a real-space survey, i.e. a survey which plots galaxies at their conformal distances. Then, 
we estimate the sensitivity of these statistics in the linear regime for a survey with depth and geometry 
comparable to the SDSS by determining those models which can be excluded to given confidence level when 
a particular fiducial model is assumed. In this way, we hope to put into perspective the ability of upcoming 
redshift surveys to fix cosmological parameters in the absence of external constraints. 
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2. Redshift Distortions in The Power Spectrum 

The theory of the hnear distortions of redshift surveys with regards to the underlying real-space 



distribution was first investigated by Kaiser (1987) for shallow surveys where redshift evolution and 



geometrical effects were negligible. This work has been extended to include both evolution and geometry 



for larger redshifts by Matsubara fc Suto (1996) and Ballinger, Peacock, fc Heavens (1996) ; for convenience 



and clarity, we shall repeat their work following the formalism of Matsubara and Suto. We begin with the 
standard assumption of a Friedman-Robertson- Walker universe with a metric given by 

ds^ = -df + a(tf{dx'^ + S{xf{de^ + sin^ 9 d./.^)}, (1) 

where S{x) is determined by the geometry of the universe through the spatial curvature K, and 

sm{^/Kx)/VK iK>0), 
S{X) = { X iK = 0), (2) 

sinh(V^x)/V^ (K < 0). 
The curvature is given in terms of the present day Hubble constant Hq, the matter energy density f2o and 
the energy density in cosmological constant Aq; 

K = H^iQa + Xo-l), (3) 

where we assume that the present scale factor oq is unity. The proper radial distance from an observer to a 
source, Xj can be determined from the integral 

x=r— (4) 

^ J, a{t) Jo H{zy ^ > 

where in the final integral, we introduce the redshift-dependcnt Hubble parameter 



H{z) = i/ov/^^o(l + zy + (1 - Oo - Ao)(l + z)2 + Ao. (5) 

To see how redshift-space maps are distorted with regards to their real-space counterparts, we need to 
consider both the geometry of the universe and the peculiar -velocities of the objects being observed. Let 
us first consider geometry. We would like to examine small displacements about a given origin located at 
a redshift z with respect to a terrestrial observer [z ^ 0). This small displacement is represented by the 
vector X in comoving, real-space coordinates and has components x\\ parallel to the line-of-sight and xj_ 
perpendicular to the line-of-sight. To first order in the Taylor expansion we may write 

x\\ = oz — —^oz 

" dz Ho 

xi_ = s{x{z))5e^^^5e, (6) 
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where C|| = Hq/H[z) and c±_ = HoS{x{z))/ z. We have also assumed a distant observer by Unearizing m d( 
Since cn and c±_ have different functional forms, a spherical object will be distorted into an ellipse when 



mapped into redshift space, and it is this effect that [Alcock fc Paczyhski (1979)| suggested could be used to 
measure cosmological parameters. 

For galaxy surveys we also need to consider how the peculiar velocities of objects affect redshift 
distributions. A galaxy that would be observed at redshift ztme in the absence of any peculiar motion, 
and which has a non-relativistic peculiar velocity v relative to the background, is seen by an observer with 
non-relativistic peculiar velocity iP with an apparent redshift 



1 + Zapp = (1 + ^t™e)(l + V\\ - V?,). 



(7) 



Here U|| is the peculiar velocity of the object projected along the line-of-sight. Using eq. (Q) along with the 
eq. we can calculate the displacement in redshift-space s for a given x including both the geometric 
and velocity effects. If we drop the term proportional to + f y — v^), which is small in comparison to 

the others, we can easily show that 



si = 



S3 



c^{z)' 



S2 = 



X2 

c^{z) 

1 



l + z 
H{z) 



Ho C||(z) 

where we recall that z is the redshift of the origin of reference and ztme — ^ is equivalent to 5z in eq. 
(^. Note that the 3 direction points along the line-of-sight. The above expressions now give us almost 
everything we need to know to relate the real-space density-contrast to the redshift-space contrast. The 
conversion is made by considering the Jacobian transformation from x to s: 



(8) 



dx 

ds 



1 - ^T^—V\\[X) 



(9) 



H{z) dx3 

to linear order in the velocity perturbation, where the local density in real space S^{vecx) is thus enhanced 
by this factor when observed in redshift space. For the mean density, we take the ensemble average, so the 
velocity term vanishes when averaging the Jacobian. Putting everything together, we get 

l + z d 



S'{s{x)) = 6''{x) 



-v\\{x). 



H{z) dx3 

This result is true for any fluid component, e.g. galaxies, dark matter, etc., as it depends only on the 
continuity of the fluid, but in the special case of the total mass fluctuations (5™ we can go further and 
replace the velocity term by assuming, in the linear regime, that fluctuations grow only by gravitational 



(10) 
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instability. Perturbation theory gives us a relationship between (5™ and the peculiar velocity, namely 



v\\{x) = -^J{z)d,A-H-^{x), 



(11) 



where A ^ is the inverse of the Laplacian operator (see e.g. Peebles 1980) and f{z) is the logarithmic 
derivative of the linear growth rate D{z), 



d\iiD{z) 
din a 



VL{zf 



A(z) 
70 



n{z) 



In this equation, the redshift dependent cosmological parameters ^l{z) and A(z) are given by 



Q.{z) 



Ho 



H{z) 



{l + zfno, A(z) = 



Ho 



H{z) 



(12) 



(13) 



(during the matter-dominated epoch) while the linear growth factor has the well known solution (Peebles 

1 + z' 



(1980)) 



D{z) = '-^H{z) 



H{z'f 



dz' . 



(14) 



Galaxies do not necessarily trace mass, but, to first order we may assume that the relationship between 
galaxy and mass fluctuations is linear, given by a bias factor 5^ — b{z)6™'. The redshift dependence of 
this bias factor can be quite complex, depending on the dynamics of galaxy formation, but for z ^ 1 it is 



reasonable to follow Fry (1996) and suppose that all the galaxies in our survey were formed well prior to 
their observation with some intrinsic bias. Treating these galaxies as a separate matter component in the 
perturbation equations. Fry has shown that 



(15) 



Putting everything together, we may now write the linear overdensity function as a function of s in the 
neighborhood of a given origin at z: 



0x3 



(16) 



where P{z) = f{z)/b{z). Both for expressing the density fluctuations, usually assumed to be a Gaussian 
random field, and for evaluating the inverse Laplacian, it is convenient to work in Fourier rather than real 
space. Thus we can write 



(2^ 



6(0)i^(0)'^ 



(17) 



where 5{k) is the Fourier transform of the real-space density contrast evaluated at z = 0. 
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The above result is more than adequate for perturbations on the largest scales, but on smaller scales, 
we must account for the non-linear evolution. This causes the density fluctuations to rapidly virialize into 
halos, and the effect on the redshift-space density field is well approximated by adding a random velocity to 
all the mass elements. Empirical evidence suggests that the distribution of these velocities is well described 
by a power-law model with dispersion (|Cole, Fisher, fc Weinberg (f995j ). Thus, the corrected density 
fluctuation in redshift space then takes the form 



(2^ 



1 + 



n -1 



b{z)D{z) ~ ^^.^(^ 



(18) 



Now let us consider the Fourier coefficients S'^{K), which we would derive from the redshift-space density 
distribution, defined as 

S'{K)= [ d^s 6'{six))e-'^ '. (f9) 



Substituting in the definition of (5*(s(x)) from eq. (|18|), we get 



_fk_ 
(2^ 



1 + Piz) 



fc2 



b{z)D{z) ^r^f„yk-x{s)^^tK-s 



mm 



Recall that x{s) = c±s± + cush, so by performing the integration over d^s we are left with 



l + P{z) 



fc2 



mm 



(20) 



(21) 



where S^{K — c±k± + C||fc||) is a Dirac delta function. The final integration is now trivial to perform, and 
doing so leaves us with 



1 



,-2 -2 



1,1 -2 2 7^2 2 
1 + 2*^11 '^^^ ^ 



X5r (ify^(c^2_ ^-2)^2+^-2^ (^n^S^) ^ 



(22) 



where fj, — K^/K. 



3. Fluctuations in Galaxy Counts 

In the previous section, we derived an expression for the Fourier components of the galaxy density 
field mapped into redshift space valid in the neighborhood of some observation point. Unfortunately, the 
results fail as s grows large enough to break the linearization constraints enforced to calculate c\\ and c^, 
so we can only apply them to small volumes. The question becomes, how can we best use these results to 
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analyze redshift surveys which cover large angles of sky and are deep in redshift space? The only practical 
alternative is to break up such a survey into many sub-volumes, and measure the distortions within each of 
them. These measurements can then be fit to theoretical calculations and used to determine cosmological 
parameters. In |3.l|, we shall discuss the calculation of multipole moments of the power spectrum as 



observed in a single sub-volume, while in §3^ we will calculate the cosmic variance of these moments, for a 



single volume. Finally, in §3.3, we calculate the shot noise or finite sampling variance for this sub-volume. 
Together, these computations will allow us to determine the sensitivity of redshift surveys which contain 
many statistically independent sub-volumes. 



3.1. Moments of the Galaxy Distribution 

The mean value for multipole moments of the power spectrum derived from a finite sample window 



were first calculated in Cole, Fisher, & Weinberg (1994); here we repeat their calculation, to which we will 



add a calculation of the variance of these statistics in the following subsections. Since we are considering 
galaxy surveys as our source of data, we must consider not continuous density fields, but discrete realizations 
of those fields given by the point locations of the galaxies. This is a well understood problem, so we shall 



follow Peebles (1980) , and subdivide our survey into a large number of cells such that the probability of 
finding two galaxies in a cell is vanishingly small in comparison to that of finding one in a cell. We define 
Ni to be the number of galaxies in cell i, where by our supposition iV^ is either one or zero. In the previous 
section, we considered the Fourier components of the density fluctuations in a neighborhood about a fixed 
redshift. For a corresponding survey, we ought to only consider galaxies inside a window such that the 
linearization performed previously is a valid approximation, and thus we will weight each of our cells by 
a window function w(si) with volume Vw Also, real galaxy surveys usually do not uniformly sample all 
the volume in the survey, £.17. a flux limited survey sees fewer and fewer objects as distance increases; the 
function (/'(s'i) will represent this selection effect. Putting everything together, we see that a sensible way to 
calculate the Fourier coefficients of the density fluctuations is 

h = - (^^)) ^tSv^"^"^-"' (23) 

nVw ^ (p(Si) 

where n is the mean number density in the volume and angle brackets refer to the ensemble average. Dividing 
by (j) will, as we shall see, cancel its effects on the statistics we calculate rendering selection-independent 
results. The expectation-value for the number of galaxies in cell i is given by (Ni) = n(f){si)(Psi. As we 
shall see momentarily, dividing the 5^ 's by the selection function will remove the effects of the latter from 
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the statistics we calculate. 



In the previous section we showed that the Fourier coefficients of the density field in redshift space are 
dependent on /i, the cosine of the angle formed by K and the observation axis. Furthermore, the angular 
dependence is a function of cosmology which appears through /3, C|| and c^_. One might hope Nakamura J 
Matsubara, & Suto (1997) that by measuring the angular dependence of the (5^'s one can determine the 



underlying cosmology, and thus, to that end, we consider the multipole moments of the square of the 
Fourier coefficients 

p,{K) ^ j dQKViifiK) {m)~6iKr) , (24) 

where Ve is a Legendre polynomial. We choose the square of the Fourier coefficient because the ensemble 
average of the coefficient itself is zero and therefore we would see no signal in the mean. Substituting in for 
Sj^ from eq. (^3|) and taking the ensemble average, we find 
2£+l 



An 



<P{Si)(j){Sj) 



The remaining expect at ion- value can be written 



where dij is the Kroniker delta, ^ij is the correlation function given by 



_fk_ 

(27r)3 



(25) 



(26) 



(27) 



and P{k) — 6^ {k)S'^* (k) / {2tt)^ is the redshift -space power-spectrum derived from the perturbation spectrum 
given in eq. It may be conveniently expressed as an expansion in multipoles 



p{k) = j2vMpm, 



(28) 



where = k^/k is the cosine of the angel formed between k and the line of sight. Putting everything 
together, wc find 

2^+1 



PdK) 



n (j){si)4){sj)d Sid Sj 



(29) 



d^k 
(2^ 



P{k) + n(j}{st)d-^ StStj 



w{Si)w{Sj) 



<t>{Si)4>{Sj 



In the limit of infinitesimal volumes, the sums over i and j reduce to integrals, and we see in the first term 
that these spatial integrals give us Fourier transforms, reducing the above equation to 



2e + i 

An 



dflKVeifJ-K) 



P{k) uM\K + k) 



i2n) 



(30) 
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where 



(31) 



To further reduce the integral in eq. (|30| ) , we shall assume that our window function is spherically symmetric 
implying w{k) = w*{k) = w{\k\). We can then write out the multipole expansion 



w'^^\k + k) 



[K + fc)j = + k) = Y, Vnim,k)Wi'-^\K, fc), 



where fJ.K,k is the cosine of the angle formed between K and k. More generally, we will define 



W'-P-"^ {K + k)= w'-P^ [K + k)w^'i^ {K + k) 



and use the multipole expansions 



[K + k) =J2'Pn{m,k)wiPHK, k). 



and 



W^P''i\K + k)=Y,Vn{t^KMWl?^'\K,k). 



(32) 



(33) 



(34) 



(35) 



As we have already mentioned (eq. (p8|)), P{k) can also be written as a multipole expansion in /i^. To make 
use of these expansions, we will apply a well known property of Legendre polynomials, namely 



dVlkVn{i^K,k)'P„i{ilk) 



An 



-'PinK)S„ 



2n+ 1 

Using this theorem along with the expansions of the window function and power spectrum, one can see 
after some algebra that 



(36) 



PiiK) 



Pe{k)W^^'^\K,k) + 



(37) 



where the second term vanishes as the number of galaxies in the survey volume grows large. Conversely, it 
is easily removed from measurements in surveys for which shot noise makes a significant contribution. 



3.2. Variance of the Moments 

Now that we have derived the multipole moments of the power spectrum for a finite volume, one can 
imagine taking many such volumes in a redshift survey and calculating pi's as a function of z. Then by 
using the theoretical results, one could fit various cosmological models to the data and determine the best 
fit. Our ability to fit the data is limited by two types of noise: shot noise and cosmic variance. Shot noise 
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arises because galaxy surveys do not include an infinite number of objects, but it scales as the inverse of the 
number of galaxies in the sample and is negligible when N ^ 1. The second source, cosmic variance, arises 
from the simple fact that we live in only a particular realization of the ensemble of possible universes. No 
amount of galaxy sampling can eliminate this noise, so it sets a theoretical upper limit to the accuracy of 
any measurements that we might make. In this section we will calculate the intrinsic cosmic variance that 
we can expect in our measurements of the moments of the power spectrum for a single window, representing 
the ideal limit achievable by any redshift survey. 

Let us begin with the usual definition for the covariance, namely aj^,{K,K') = 
(pc {K')pi{K)) — {p£{K)) (pc {K')) . We choose different magnitudes for K and K' so as to treat 
the most general case. Using our definition of pe we can write 

'2e + i2i' + i f ,^ ,^ ^, , 1 



<jI,,{K,K') 



Att 47r 



■ / dnKdnK'ViiiK)V{i^K')^77^ (38) 



^ ^«(■?0w(gJ^)w(4)^^(g^) ^iX■(j^.-r,)^»K^(j•,-r,) 



X ((AT, - (Ar.))(7V. - {Ni)){Nk - {Nk)){Ni - {Ni))) 

-mmpHK')). 

In the limit that nVyj grows large, i.e when ignoring shot noise, it can easily be shown that the ensemble 
average of the product of A'''s is 

n'^d^Si d^Sj d^Sk d^si (l>{si)(l>{sj)(l>{sk)(l>{si) [iijiki + iikiji + ^.iiikj] , (39) 

if the underlying distribution is Gaussian. Restricting ourselves to Gaussian fields is reasonable on large 
scales, particularly for inflationary models of structure formation. Putting things together, we get 

"2^+12f + 1 



^2 



[ d^Kd^K'V{^lK)V{^lK')^ (40) 



X [^ij^ki + ^ik^ji + ^ii^kj] d^Sid^Sjd^Skd^si 
-{Pe{K)){Pr{K')) 

The three combinations of two-point correlation functions ^ij^kl, ^ik^jl, ^il^kj, give us three different terms 
to calculate, which we call /i, I2, and Is respectively. Close inspection of Ii reveals that it is exactly 
equivalent to (pe{K)) {pg/ {K')) canceling the final term. Thus we are left to calculate I2 and 73. 
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Now let us take a closer look at l2- As before, we can do the spatial integrals and replace the window 
functions with their Fourier transforms. The resulting equation is 

2£+12f + 1 f .„ .„ , f d^k (Pk' 

J2 = 

47r 4tt 

xw^^\K + k)w^^\K - k')w^^\K' + k')w^^\K' - k). 



/f rfik d'^k' 



We now expand the window function in a multipole series: 



n . 



n m— — n 



where Ynm is the usual spherical harmonic and the last line was derived from the addition theorem of 
spherical harmonics. Expanding similarly the other window functions, we may rewrite eq. as 

T = V V V V (4^)^(2^+ l)(2f + l) 

2 2^ 2^ 2^ 2^ (2ni + l)(2n2 + l)(2n3 + l)(2n4 + l) 

ni.nii 712, m2 '^z^'^^ 714,7714 



(27r)3 



(27r)3 

X / dQ.KVi{jJ.K)yni,'mi{^K)y7i2,iyi2{^K) 



X y dQ.K'Vtj{^lK')yn3,m3{^K')yni,mi{^K')- 

Note that the factors (—1)"^ and (—1)"-' arise from the fact that 7'„(— x) = (— l)"7^„(a;). From the above 
equation, we see two integrals of the form 



^ dnKVl{nK)Yn,,mA^K)Yn2,mA^K). (43) 

which evaluate to 

\j 2£ + |^^"ii^'"2;'ni'^m2,-mi I (44) 

where Cni/,n2;mi6m2,-mi — uulcss thc triangle condition \ni — £\ < ni < ni + £ is satisfied, in which case 



1 ni 


£ 


:)( 


\ 








C'ni/,n2;mi<5m2,-mi — \/ | • (45) 

477 I n n o ; \ -mi 

The arrays in parentheses are the Wigner 3j symbols. For a complete discussion of these symbols see, for 
example, Edwards (1957)| . Also, a fast recursive algorithm for calculating the 3j symbols is described in 
Schuhen fc Gordon (1975] . 
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The second type of integral which appears in eq. ( |4^ ) has the form 



(27r)3 

Expanding the power spectrum into its multipole moments, we can write this integral as 



i 

where we have defined the function 

(fc, K)wi'] {k, K')P, (fc). (48) 

Now we can combine all of the terms together to write 

^2 — / ^ / ^ ■ ■ ■ / ^ *^mi , — 7712 , — m3 ^7714 , — 7712 ^7714, — ms V / V "^^ / 

i,j Til ,"^-1 7l4 ,7774 

(4^)4 5„,,„3;,(^C, K')S„,^„,.j{K, K')V2lTlV2FTT 



(2ni + l)(2n2 + l)(2n3 + l)(2n4 + l)V2i+TVW+^ 

^ C'm ,2,n3 ; — mi ^'7^2 j',n4 ; — mi C'yj^ ,£,n2 ;mi C'^^^ ,714 ;m3 ; 

with the implicit constraints 

\ni-e\ <n2< ni+£, 

\ni - nsl <i< ni+ n^, 

\n3-£'\ <n4< na+f, 

1^2 - n4| < j < ^2 + Hi. 

Applying all of the constraints and summing out the Kroniker delta's, I2 takes the final form 

ni+l ns+i' ni+ns n2+»4 "l 

= E E E E E E (50) 

ni.n^ n2 — |77i— ^1 774^1713— i^|7ii— 713I j — |7i2— ^^4| 777i — — 7ii 



(2ni + l)(2n2 + l)(2ri3 + l)(2ri4 + l)\/2TTTy2jTT 

^ ^ni ,i,n3 ; — mi C'^^g ,j.n4;— rni ,^ ,712 ;mi ^'7^3 ,£',714 ;mi ■ 

A similar calculation can be performed for for I^, but it should be apparent from inspection that the result 
will be the same as for I2 except that there will no longer be the term (— 1)"^+"-'. Although calculating the 
infinite sum seems daunting, we found that, for the examples calculated in this paper, only a few moments 
in the power spectrum and a few tens of moments in the window function were necessary to get very good 
convergence, making the problem numerically tractable. 
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3.3. Shot Noise 

In the previous section we considered the variance in the muhipole moments for a single measurement 
in the continuum limit. For magnitude limited surveys, the number of galaxies contained in a given volume 
will fall with rcdshift until the variance arising from finite sampling, i.e. shot noise, can become comparable 
and then dominate the cosmic variance. To accurately determine the sensitivity of real surveys, one must 
accurately model the shot noise as well. In this sub-section, we shall discuss the shot noise calculation, 
leaving the details to the appendix. 

Recall in eq. (^8|) that to calculate the covariance cr| ^/ we needed to find the ensemble average of a four 
point moment 

{{N, - {m)iNj - {NMNk - {N^ym - im)), (5i) 

which was summed over all indices. We claimed that in the limit of large n this average could be reduced to 
eq. (^9|), and we will now show this explicitly while calculating all other terms. We can rewrite the above 
average to explicitly include terms in which one or more of the indices are equal. The result is 

m - imiNj - {NMNk - {Nk)){Ni - {Ni))) + (52) 
+ (S^jiiN, - {N,))\Nk - (Nk))iNi - (iV,») + 5 permutations) 
+ {S^J6kl{{Ni - {Ni)f{Nk - (A^fe>)^> + 2 permutations) 
+ {S^JSjk{{N^ - {Ni))^{Ni - {Ni))) + 3 permutations) 

where different indices are presumed to be unequal. Since in the limit of infinitesimal volumes 
(A^f ) = n(j){si)d'^Si, the ensemble averages reduce to 

n^(l){si)(j){sj)(l){sk)(l){si)d^Sicfisj(fsk(fsi [^yCfei + ^ik^ji + Cu^jk] (53) 
+ (5ijn'^(j){sj)(t){sk)<j){si)d^SjSskS.siS,ki + 5 permutations) 
+ {5ij5kin'^(j){si)(l){sk)d^Sid^Sk{l + 6fc) + 2 permutations) 
+ {SijSjk(f>isk)(t){si)d^Skd^siiki + 3 permutations) 
+ 5ij5jk5kin(f){sk)d^Si, 

for Gaussian fields. The highest order term in n is the one which we have previously calculated, while the 
others are small in comparison when n is large. What remains is to substitute eq. (|5^) into eq. ( |3^ ) and 
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evaluate all the terms. These details are left to the appendix, but here we give the results. To all orders in 
n, the shot noise contribution to the variance of the multipole moments of the power spectrum is 



,{K,K') = v/(2^+l)(2f + 1)1 E E 



'shot £.£ 



(47r) 



5/2 



ij.nij ni = \e-i\ n2 = \l'-j\ 



2 (— 1)-'+™^] 



(54) 



\/2r+T(2ni + l)(2n2 + l)(2j + l) ^ ' ' ' 

^Tii ,n2 ;iCni ,i.n2:mj Cm ,i,i; — mj Cn2 ,n2;Tnj 

+ S^oS.o^ I 0T^(2.2)(fc)P„(fc)+5«, [l + (-l)^] -l^w(^.^)^K,K') 



+ 



V 



(2f+l)(2rH-l,X:/^P„(.,/f!^GW(.-'^ 

2 



X ji{Ks)jt(K's)jr,{ks) 
2 47r 




+ i 



l-l'+n 



+ <5, 



+ 



iV2 2i' + 1 
2 Ait 



N^2e+1 J (27r) 



^H^f^)(i.',fc)PHfc) 
^'*.W^f^)(i^,fc)P,(fc) 



where 



(55) 



and w'^P^k), M/(P'«)(fc) and VK/'''''^(fc) are defined above. In calculating these terms, we have assumed that 
the window function is spherically symmetric (implying G{s) ~ G{s)), and also that the selection function 
(j) is constant over the window, leading us to define N — nc/fV^. 

Since factors of n now appear explicitly in the shot noise terms, we require an estimate of the number 
density of galaxies as a function of redshift. The SDSS is an r-band magnitude- limited survey, so we need an 
appropriate galaxy luminosity function. An analysis of the same band for the galaxies in the Las Campanas 



redshift survey performed by Lin et. al. (1996) has shown that the galaxy luminosity function is well fit by 
a Schechter function with M* = -20.29 ± 0.02 + 5 logh, a = -0.70 ± 0.05, and 0* = 0.019 ± 0.001 . 
For our models we will presume that these parameters will also well describe the SDSS and take h — 0.7, 
consistent with the fiducial models we shall consider. To get 10® galaxies in the total survey volume, we 
choose an apparent magnitude limit of m = 17.6 consistent with the limit rn sa 18 quoted for the SDSS. A 
second map has also been proposed for the SDSS data that would generate a volume limited survey of 10^ 
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bright red galaxies with an approximate depth oi z = 0.5 (Szalay (1997)). We found in our calculations that 
the density of galaxies for this survey was insufficient to significantly strengthen the constraints derived 
from the proposed magnitude limited survey, so we shall only analyze the latter in detail. 



4. Sensitivity of Redshift Surveys 

The next generation of redshift surveys, particularly the Sloan Digital Sky Survey, will probe large 
fractions of the sky observing up to a million galaxy redshifts. In this section, we would like to estimate 
how well SDSS, the largest proposed survey, can measure cosmological parameters from redshift space 
distortions. Specifically, we want to determine whether the geometric and evolutionary effects which occur 
at larger redshifts can break the degeneracy between the matter density and the bias, producing a clean 
signal for f2o- To that end, we consider a pair of fiducial models which have good concordance with 
observation and perform a statistical fit of cosmological parameters. 

In the last section, we calculated the mean and variance of the multipole moments of the power 
spectrum measured in small sub-volumes of a large survey. For nearby samples, when geometric effects 
are unimportant, the ratios of any two moments depend only on f3 and the shape of the power spectrum, 
but not its normalization. Cosmological models are roughly degenerate in parameter choices for which 
/3 — r^Q-^/feo is a constant. We shall reconsider these ratios with the inclusion of geometric and evolutionary 
effects to see if this degeneracy can be broken in deeper surveys. The ratio of P2/P0 is considered as it yields 
the largest signal to noise proportion; our numerical results show that the non-linear clustering effects tend 
to suppress the signal in higher moments. Two types of fiducial models with fig = 0.3 are considered, 
one an open model with A = and the other a flat model with A = 0.7, as these choices are favored by 



observations of large scale structure, e.g. Tadros & Efstathiou (1995). For the dark matter component, we 



assume Cold Dark Matter, with a fixed shape parameter of F = ilo^o = 0.2, which we use for all test models 
so as to deconvolve the redshift distortion effects from the effects of the of the changing power spectrum 
shape. Furthermore, F should be well measured in the upcoming redshift surveys by direct calculation of 
the power spectrum, and thus should be a fixed input. To model the non-linear dispersion, we choose a 
velocity of 300 km/s which corresponds to a length scale of 3 Mpc in Hubble units. 

In comparing our test models to the fiducial ones, we must consider how one can best extract 
measurements of p2 /pa from a real survey. Cole, Fisher, fc Weinberg (1995)] examined n-body simulations 



in the case of no evolution or geometric effects. They took repeated sub-samples by randomly locating the 
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center of a window function and measuring the multipole moments in these sub-volumes. To ensure that 
the data are statistically independent, the typical separation distance between sub-samples is chosen to be 
the size of the window, otherwise, the overlap will introduce correlations. We reconsider this procedure, 
only now, the effects of redshift evolution and geometry will be included to more accurately model real 
surveys. 

Clearly, to maximize the number of measurements, one ought to consider the smallest possible 
sub-volume. Unfortunately, non-linear effects dominate on small scales and our simple model of the 
non-linear velocity dispersion will break down. In n-body simulations, [Cole, Fisher, fc Weinberg (1995) 
found that the corrected linear models failed at wavelengths below about 20/i~^ Mpc, so we shall take 
this to be the lower limit on which we can apply our linear calculations. Thus we consider a Gaussian 
window — chosen because its multipole moments are expressible analytically — with a radius in redshift 
space which corresponds locally to tq = 20h^^ Mpc, where r = Hq^z. For our calculations of pi[K), we 
selected K = 27r/20 /iMpc^^ consistent with the n-body results. We then divided space into slices -Kz^dz 
corresponding to the tt steradians to be covered by the SDSS, with dz equal to the window diameter of 
40/i~^ Mpc locally. Redshifts corresponding to less than twice the window width were ignored because 
they were not sufficiently distant for the distant observer (i.e. small angle) approximation assumed in eq. 
(^ of §H to be valid. High redshift data was dropped once it became shot-noise dominated and no longer 
contributed significantly to the fits. In each of these slices we divided the shell volume by the window 
volume to determine the number of statistically independent measures that were available at that redshift. 
To determine the expected variance of the p^{K)'ii, we divided the variance for a single measurements by 
the number of volumes, the standard suppression for multiple independent measurements, and calculated 
the errors in their ratios using simple propagation of errors. In other words, the error bars were determined 



where M is the number of independent volumes in the particular redshift slice. The result was a set of data 
points showing the expected mean and deviation for the redshift bins one would reasonably choose when 
analyzing a real survey. We then repeated this procedure for tq = AQh^^ Mpc, the smallest scale which is 
independent of the 20h^^ Mpc data. We found that data on both scales were necessary to constrain all the 
model parameters, while including larger scales was not significantly more constraining. 

To compare with other cosmologies, we calculated the mean values of P2/P0 for models with varying 
rio, 60 1 and a and calculated a using the fiducial model. We defined our confidence limit as the surface of 



by 




(56) 



- 18 - 



constant into which a given percentage of best fit model parameters would fall for multiple realizations 
of the fiducial ensemble. For example, the 95% confidence limit is defined such that the best fit model 
parameters will fall within that surface 95% of the time. Since we have three fitting parameters, there are 
three degrees of freedom in our statistics. 



5. Results and Conclusions 

In figure ^ we show two examples of P2/P0 plotted as a function of redshift, where we have removed 
the shot noise contributions to the mean values. In the upper panel, we show a fiducial model (stars) 
of ilo = 0.3, Ao = 0.7 and 60 ~ 1-0 with la error bars, compared to a test model (boxes) of fio = 0.48, 
Aq = 0.52, and 60 = 1-2; and in the lower panel, we show a fiducial model of VIq = 0.3, Aq — 0, and 60 = 1.0 
compared to a test model of = 0.85, Aq = 0, and 60 = 1.60. In each example, data for both the 2{)h^^ 
Mpc and AOh^^ Mpc windows are shown, where the larger window produces higher values for the ratio, and 
all models use the fiducial value for a. Data was terminated at the redshift after which shot noise dominated 
the result. Both comparison models are barely accepted at the 95% confidence limit when the fiducial 
velocity dispersion was used for each. We emphasize this with figure ^ which shows the 68%, 95%, and 99% 
confidence limits for flat, cosmological constant test models given a fiducial VIq — 0.3, Aq — 0.7, and 60 = 1-0 
model like that of figure |l|. These surfaces are projections of the full three dimensional confidence volume, 
and they represent the confidence limits when the non-linear velocity dispersion is unconstrained by other 
data. The dotted line shows the curve rip'^/feo = 0.49 which is the degeneracy one expects locally when 
redshift effects are ignored. We do not show a figure comparing open models to our fiducial cosmological 
constant model, as none were acceptable at even the 99% confidence limit. In figure |^ we switch to an open 
f2o = 0.3, Ao — 0, and &o — 1-0 fiducial model and test it against other open models. Finally, in figure^ we 
test this fiducial model against fiat cosmological constant models. 

In figure]^ we see the significant redshift dependence of P2/PO1 where we would expect none in models 
which ignore geometry and evolution. Naively fitting a horizontal line representing models with no redshift 
dependent effects, we get a best fit (3q — 0.45 for our fiducial flat model and /3o = 0.26 for our flducial open 
model using only the 20/i~^ Mpc data; we expect a value of /Jq ~ 0.5 in both cases. Thus we verify that there 



will be systematic errors in models which ignore redshift evolution effects, a point emphasized by Nakamura, 



Matsubara, fc Suto (1997)| . The slopes of theses data can be qualitatively understood by considering the 



competing effects of the evolution of (3 and geometry. At deeper redshifts, ri(z) 1 and /3 increases. 
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tending to increase the value of the ratio of p2 /po for the power spectrum and non-linear dispersion we 
have chosen. Geometry, on the other hand, pushes the effective scale which is probed to smaller scales and 
thus larger values of the wave number, i.e. objects appear more spread out at high rcdshift. Non-linear 
dispersion tends to drive the multipole ratio to lower values as the effective K grows, so geometric effects 
tend to decrease the multipole ratios. Thus we have two competing effects which determine the slope of the 
multipole ratio curve as a function of redshift. For Qq — 1 models, /3 remains fairly constant — changing only 
due to the evolution in bias. These models are dominated by geometric effects and thus have the steepest 
negative slope. The high cosmological constant model has the flattest slope, because the rapid evolution 
in f3 just cancels the geometric effects, while the open models fall in between. The effects of changing the 
non-linear velocity dispersion generally shifts the overall normalization of the data without affecting the 
slope significantly. With only the 20h~^ Mpc data, there is a degeneracy between bo and a. To break it, 
we need to include the data from the larger window, as the change in normalization is different for the two 
window scales. We also note that the nearly flat evolution of the multipole ratio also explains why open 
models cannot fit our fiducial cosmological model; all are too sloped to be good fits to the data. 

Regarding the determination of cosmological parameters, the results shown in the last three figures 
reveal interesting results. For the fiducial open model (figure the degeneracy between models of differing 
r^o and bo is weakly broken, and the redshift distortions do permit us to distinguish between the open 
fiducial model and fio = 1 models at the 95% limit, although not at the 99% limit. However, when we 
tried a fiducial ilo = 0.4 model, = 1 was allowed at the 95% limit, so one can conclude that an open 
universe with Qq = 0.3 is just on the edge of being able to exclude critical matter density models. For the 
case of the flat fiducial model (figure ||), we see that values of flo > 0.48 are ruled out at the 95% confidence 
limit, producing significantly stronger constraints than is the case for the open model. In figure ^we test 
the likelihood of confusing flat models with open. We see that models with some cosmological constant 
can be confused with the open fiducial model. Overall, we conclude that redshift surveys like SDSS may 
just be able to determine bias and evolution independent measures of cosmological parameters, at least in 
discriminating the extremes of fig = 1 and Qq = 0.3. although only nominally so for open universes. We 
also see that naive estimates of /? which ignore redshift evolution are systematically biased towards smaller 
values. 

When it comes to determining cosmological parameters, the ideal result is to measure them with 
several independent observations, looking for a consistency of result which would demonstrate that we 
understand the fundamentals of cosmology or an inconsistency that would indicate a failure of some aspect 
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of current theory. It has been suggested that redshift surveys offer promise of being a rehable source 
of data for such cosmological parameter determination, when appropriate statistics are apphed. Simply 
measuring the change in the mean number of galaxies as a function of redshift is inadequate, because we 
cannot observe the gravitational component in the absence of reliable models for the non-gravitational 
evolution. Alternatively, one can consider the redshift evolution of the multipole moments of the linear 
power spectrum, which are driven only by gravity, and hope to cleanly observe cosmological effects. Our 
results suggest that if the true model of the universe contains a large cosmological constant, then Qq is 
tightly constrained; however, for open models the the extremes of ilg = 1 and Qq = 0.3 only can marginally 
be distinguished (95% but not 99% confidence limit) in surveys on the scale of SDSS. 
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Appendix A. 



In §3.3, we discussed the complete expansion of the variance including the lower order terms in n. Here we 
shall show the details of calculating the shot noise terms ignored in § |3.2| . We begin with the leading order 
correction in eq. (p2[), which when substituted into eq. (pq), may be written 



2/ -I- 1 2/' + 1 /■ 1 



X [%((iV, - imfiNk - {Nk)){Ni - {Ni))) + 6,k{{N, - {N,)f{N, - {N,)){Ni - (Ni))) 
+ Su{{N, ~ {N,))'{N, - {N,)){Nk - (iVfe))) + S,k{{N, - (A^.))(iV, - {N,))'{Ni - {Ni))) 

+ s,im - (mm, - {N,))'{Nk - (iv^))) + Skim., - {mm, - {N^mm - {m))')] 



Now let us digest this term by term. The first we shall designate Ii i and, recalling eq. (p3D, we may write 
it as 



X 



To allow us to evaluate many of the integrals that we will encounter in this section, we need to make a 
simplifying assumption about the selection function 0: that it is approximately constant in the window 
of interest. For small windows this should be reasonable, so we thus define N ~ n(f>Vw to be the number 
of galaxies in a given window. Substituting in the Fourier expansion for the correlation function and 
evaluating the spatial integrals, it is straight forward to show 



= Sio-— I dnK,ViifiK') I 77^W^'''\k + k')P{k), (3) 



with 

V ^ J d'^sw^is). (4) 

and defined in eq. (|3^). Expanding everything into Legendre series and evaluating the angular 

integrals, one can show directly that 

4n V f k^dk 



(5) 
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By inspection of eq. (|l|), one can also see that Ii^ is equal to Ii^i under the interchange of £ and £'. However, 
looking at eq. (^7|), we see that both of these terms will cancel with the terms coming from {pi){pt), so 
they are dropped below. 

Now we continue on to consider the second term I12 which may be written 

Performing the spatial integrals, we find 

1 2/' + 1 2/" + 1 /■ 

h.2 = J— / dnKdnK:V,{llK)Vi>{tlK') (7) 

w^^^ {K + K')w{k + K)w{k - K')P{k). 



(27r)3 

with w^^) defined in eq. (^l]). Let us take a closer look at the integral 

d^k 



,^ P{k)w{k + K)w{k^ K'). (8) 

If we expand each term into a Legendre series, and then expand those Legendre polynomials into spherical 
harmonics using the addition theorem, we can see after evaluating the angular integrals that this term may 
be reduced to 



ni+n2 ni /. n2 



(2ni + l)(2n2 + l)'^"^'"^'*^ ^Y2j^J^-^"i>— '"i(^-^-'-^"2,mi(f^if')C'ni,i,n2;mi (9) 



('2^1 4 

ni,"2 i=|ni— n2| mi= — ni 

with 5 defined in eq. (|4|). To make further progress, we substitute the above back into eq. (^) and expand 
if)^^-' {K + K') in a spherical harmonic series. The rest of the work consists in evaluating various angular 
integral like we have already seen, so we shall just quote the final result: 

= V(2^+1)(2£' + I)l E J2 E (-1)"^+"^- (10) 



(4^) 



N 

5/2 



V2r+T(2ni + l)(2n2 + l)(2j + 1) 

X Wj (-^; jt^ni ,n2 ,i,n2 ;mj ^ni ,£,2; — mj ^n2 ,n2 ;mj ■ 

By inspection one can see that the term Ji_5 is equal to The remaining terms /i,3 and Ji^ are 

equivalent to /i_2 under the transformation — 1"^ — > —1-'. 
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The next correction term in eq. (|5^) has the form 

2^ + 1 2^' + 1 



An 47r 



dnKdQK''P{p.K)'P{p.K' 



(11) 



w{si)w{sj)w{sk)w{si) 



»ffc,i <t>{si)(l>{sj)(l){su)(t>{si) 



(12) 



We refer to each of the terms here as l2,i, ^2,2 and 12,3 respectively. The first, /2,i can be expanded into 
the following: 

+ /(2^)3 

after we insert the Fourier expansion for the correlation function and evaluate the spatial integrals. The 
first term will cancel with pieces from {pe){pe'), leaving the second. The angular part of the k integral is 
easily evaluated, leaving 

h.i = SeoSi^oj;^ j ^W^d'D (fc)Fo(fc). (13) 

For the second term in eq. ([Ill), we shall not automatically evaluate the spatial integral because the 
results are simpler to reduce for some terms if we evaluate the angular parts first. Thus we write 

1 9^ 4- 1 op' 4- 1 r . . 

dnKdn'^ViiiiK)vi{p'K)Y. 

•' k,l 



iV2 Att Att 



Si 



(14) 



The first term, designated I2.2A evaluates to 

1 2^ + 1 2£' + 1 



1 



d^k 
(2^ 



I2,2A — 



iV2 Att Att 



P(k) 



drtKdnK'riitxK)ri'{^iK')w^^-^\K + K'), 



d^Skd^si. 



(15) 



which, after expanding the widow function into a Legendre series and evaluating the angular integrals, 
reduces to 



(16) 



To tackle the second term in eq. (|l4) , we first define a new set of spatial variable s = Sfc — si and 



s' = Sk — si. We rewrite the spatial sums as integrals in the new variables, producing 



'2,2B 



1 2i+12£' + 1 
iV2 Att Att 
1 fd^sd^ 



dVlKd9.K'Vf,{^.K)Vt'{^.K') 



(17) 



^ 2 



^-^(^)-^(^)e^^-V^ 



e'^-'P{k). 
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The quantity 



1 f cPsd^s' n.s + s', 9/S — s*. 



(18) 



is a function only of the magnitude of s if w is spherically symmetric. The exponential has the following 
expansion in Legendre polynomials 



(19) 



which we may use, along with eq. ( |36D and eq. (45) to evaluate the angular integrals. The final result can 
be written 

2 



2.2B — 



'n 



(2.+ l)(2.' + l)E/^P4^)/^G(.) 



(20) 



£ i' n 




■i+t+n 



u{Ks)n\K's)3r,{ks). 



Inspection of the third term in eq. ( p^ reveals that /2.3A = (—1)^-^2.271 while /2.3B = i I2.2B- 

Moving along, we see that the fourth term in eq. ( ^2|) is 
2^ + 1 2^' + 1 



47r Att 

E 



dnKdnK'V{^iK)'P{fiK' 



(21) 



w{s,)w{Sj)w{sk)w{si) _ ) ,K' . (^-, ) 



4){Si)4){sj)(t){sk)(t>{si) 



+ S,kSki{{N,N,) - {N,){N,)) + S,M{N,N,) ~ {N,){N,))] 



Defining ' according to eq. (35), it is straightforward to show that 

14^ f l^,,A3..i),^, 



(22) 



iV2 2^+ 1 J (27r)3 

The term 2 — -^3,1 while the terms 13^3 and 13^4 are equivalent to i under the transformation i £' and 
K ^ K' . The final remaining term in eq. (|38| ) can be evaluated trivially, producing 



(23) 



Putting all of these (eq. 



13 Id 2 



'Lz_il 'l~jil'rL~I1 



l§ yields eq. (|54|) 
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Figure Captions 

Fig. 1. — Comparing the redsliift evolution of P2/P0 for different cosmological models. In the upper panel, 
the starred data points show our fiducial flat model (Aq = 0.7, 60 = 1) with error estimates compared with 
the square data points which represent a Aq = 0.48, 60 = 1-2 model. Data for windows of 20/i~^ Mpc (lower 
curves) and 40/i~^ Mpc (upper curves) are shown. The comparison model has the smallest value of Ao which 
is still acceptable at the 95% confidence limit. In the lower panel, the starred data represents our fiducial 
open model (fio = 0.3) compared with an = 0.85, 60 = 1-6 model, which is again just acceptable at the 
95% confidence limit. 

Fig. 2. — The 68%, 95%, and 99% confidence limits for l^o and 60 with a unconstrained, when comparing 
flat, cosmological constant models with a flducial Oq = 0.3, Aq = 0.7, a = 300 km/s, and 60 = 1 model. The 
dashed curve plots ClQ^/bo = .5, the naive degeneracy expected 

Fig. 3. The 68%, 95%, and 99% confidence limits for Qq and bo with a unconstrained, when comparing 
open models with a fiducial Oo = 0.3, Aq = 0, cr = 300 km/s, and bo = 1 model. The dashed curve plots 
flo^/bo = .5, the naive degeneracy expected 

Fig. 4. — The 68%, 95%, and 99% confidence limits for fio and bo with a unconstrained, when comparing 
flat, cosmological constant models with a fiducial Oq = 0.3, Aq = 0, ct = 300 km/s, and 60 = 1 model. The 
dashed curve plots ^0^/60 = -5, the naive degeneracy expected 



